Multiband Electronic Raman Scattering in Bilayer Superconductors 
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A theory of electronic Raman scattering in the presence of several energy bands crossing the 
Fermi surface is developed. The contributions to the light scattering cross section are calculated 
for each band and it is shown that the cross section can be written in terms of the sum of the 
single band contributions and a mixing term which only contributes to the fully symmetric channels 
{A\g). Particular emphasis is placed on screening in bilayer superconductors. Since any charge 
fluctuation with long range character in real space is screened by the Coulomb interaction, the 
relevant fluctuations in a single layer case are induced between different parts of the Fermi surface. 
In a single band d-wave superconductor the scattering at energy transfer twice the maximum gap 
Amaa; is dominated by those parts of the Fermi surface where A^(k) is largest. As a consequence, the 
fully symmetric (v4ig) scattering is screened. In the case of a bilayer superconductor however, the 
charge transfer is possible between layers inside the unit cell. Therefore a formalism is considered 
which is valid for general band structure, superconducting energy gaps, and inter-layer hopping 
matrix elements. The spectra is calculated for La 2:1:4 and Y 1:2:3 as representative single and 
bilayer superconductors. The mixing term is found to be negligible and thus the response is well 
approximated by the sum of the contributions from the individual bands. The theory imposes strong 
constraints on both the magnitude and symmetry of the energy gap for the bilayer cuprates, and 
indicates that a nearly identical energy gap of dj.2_y2 symmetry provides a best fit to the data. 
However, the A\g part of the spectrum depends sensitively on many parameters. 

PACS numbers: 74.20.Mn, 74.60. Ge, 64.60.Ak, 74.40. +k 



I. INTRODUCTION 

In the recent few years Raman scattering has been 
proved to be a very powerful tool to study the 
anisotropy of the-superconducting energy gap in high-Tc 
superconductorsETcl. Even if Raman scattering (in clean 
materials) can give information only about the absolute 
value of the superconducting gap, the spectra clearly in- 
dicate the extremal points of the gap as the maxima and 
minima which show up as singularities in the electronic 
density of states. The spectra show that the electronic 
part has essential density of states below the maxima as- 
sociated with the maximum of the energy gap. On the 
other hand, the data clearly shows that there is only one 
sharp maximum in the spectra, indicating that the ab- 
solute value 0|£ the maximum and minimum of the gap 
are very closed" 13. Thus the gap smoothly interpolates 
between A^aa: and — Amaa and no other singular part 
exists. Such a gap occurs in d— wave superconductors. 
One can clearly state that if the gap contained mixtures 
of different symmetries, the d— wave part at least must 
be the dominant ona3. 

The experiments carried out with different polariza- 
tion orientations pick up the contributions to the light 
scattering on different parts of the Fermi surface. For 
example, incident and scattered light polarizations e^''^ 
aligned along x + ij , x -\- y or x + y, x ^ y, select then the 



Aig + B2g or Big part of the spectra is measured, where 
the X, y coordinate system is locked to the Cu02 planes. 
Measuring the Big spectra provides information mainly 
concerning the light scattering in the neighborhood of 
the kx and ky axes, B2g spectra probes mainly along the 
diagonals, and Aig is a weighted average over the entire 
Brillouin zone. The sharpest spectra are observed in the 
Big channel and thus the maxima of the absolute value of 
the energy gap must lie along the axes, supporting a gap 
of c?2,2_y2 symmetry. The recent experiments carried out 
in different hole doped high-T^ materials provide almost 
identical spectra showing very characteristio^ifferences 
between the various polarization orientationgj. 

In forming the Raman spectra, long-range Coulomb 
screening plays a very important role. Namely, those 
parts of the spectra which are coupled to the Coulomb 
fields are screened and do not appear in the limit of small 
momentum transfer (7 0, as is always the case. The 
long-range Coulomb forces play an important role if there 
is charge transfer between unit cells at large distances. 
Thus only the charges produced by charge transfers in- 
side the unit cell determine the measurable spectra. This 
transfer is between the different atoms in the cell or the 
redistribution of the electrons between different parts of 
the Fermi surface or between different Fermi surfaces. In 
this way, only the Aig part of the spectra can be coupled 
to the Coulomb forces. This argument is based on the as- 
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sumption that only one Fermi surface is relevant, which 
is very reasonable for materials with a single CuO»,plane 
in the unit cell. Recently Krantz and CardonaP have 
raised the relevant question how a double-sheeted Fermi 
surface, as occurs in materials with more than one Cu02 
planes in the unit cell, changes the above argument. 

The general idea of charge transfers between different 
regions of k space with energies nearby the Fermi level 
is very similar to the case of semiconductors, where the 
electrons or holes can be redistributed among several val- 
leys of the conduction bands or the different parts of the 
valence bandda. 

Except for scattering in the Aig channel, the presence 
of the double layer does not affect essentially the theory 
of Raman scattering as the Coulomb screening does not 
play a role. In the Aig symmetry, hflwever, additional ef- 
fects may occur as expected by RefJj. Therefore, the Aig 
channel will be examined in great detail in the present 
paper. 

The main question is whether or not in the double 
layer compound a singularity can occur at 2Ajnax in the 
Aig symmetry, similar to the one which occurs in the 
Big case. In the single layer compounds the calculations 
show that the peak with Aig symmetry is far from being 
singulaxj, with a position ranging from somewhat near 
A,„q2;E|0 to slightly below 2A„iaxa. 

The consideration will be focused on the energy range 
Lu ~ 2A,nax, thus nearby the singular contribution. If 
there is a singularity at 2Amax, then the region nearby 
the axes dominates the spectrum and from the point of 
view of screening the other parts can be ignored as their 
relative weight disappears as ui ^ 2Amax- In Fig. 1, 
those regions of are shown, which is also where the ab- 
solute value of the gap is maximal. There, the Fermi 
surface of the two layers are shown. Of course, in a more 
appropriate discussion followed in the present paper, the 
even and odd (bonding and anti-bonding) combinations 
of the wave functions of the layers must be introduced. In 
order to first get an insight the labels of the Cu02 layers 
is used. Considering the electronic density, the singular 
regions of the Fermi surface are indicated by small ellip- 
soids with a strong resemblance to the different valleys in 
the case of multivalued semiconductors. The two layers 
double the number of these regions. In the case of Big 
symmetry the charge transfer is between the region of 
the kx and ky axes, respectively, and changes sign going 
from one region to the other. In contrast, for Aig sym- 
metry, each region has the same phase on a particular 
Fermi surface. Therefore, in the case of the single layer 
material, considering the singular part, the charge trans- 
fer is only between different cells, which is screened. In 
the double layer case, the Aig symmetry allows a charge 
transfer between the two layers of both the same and 
opposite sign on each Fermi surface. Using the bonding 
and anti-bonding wave functions of the two layers, the 
odd combination plays an important role. 

Considering the interlayer charge transfer process, the 
question must be raised how the light coupled to both 



planes in a similar way can result in a charge transfer 
between the symmetrically displaced layers. That dis- 
crepancy can be resolved only by considering the cou- 
pling between the planes which automatically leads to 
the introduction of the a bonding and anti-bonding wave 
functions. The coupling leads to a split of the two bands 
which is also a measure of the possible charge transfer 
between the planes, as will be discussed in the present 
paper. A nearly singular peak at w ~ 2Amax can occur if 
the splitting of the superconducting gaps for the electrons 
with even and odd wave functions is small. The strength 
of the peak is proportional to the square of the differ- 
ence between the Raman efficiency (7+ — 7-)^, where 
7± correspond to the coupling of light to the bonding 
and anti-bonding electrons, respectively. This difference 
can be demonstrated in the limit where the incident and 
scattered photon energies are small compared to the rel- 
evant band energies. In this limit, — 7_ is propor- 
tional to the splitting of the effective mass in the two 
bands. It turns out that considering a realistic descrip- 
tion of YBa2Cu307 and Bi2Sr2CaCu208, the sphtting 
of the mass is less than 10 percentQ. Thus the peak at 
LU — 2Amax Can not have a larger relative intensity then 
10^^. If that argument holds, the Aig spectra must be 
very similar to those calculated from the two separate 
uncoupled layers and the effect of the coupling is very 
likely to be beyond the observable range. 

In our paper, we consider the most representative of 
the single and double layer compounds La 2:1:4 and Y 
1:2:3, respectively, while we will only briefly mention the 
multilayer compounds, such as Tl 2:2:3:2, in the Ap- 
pendix since the main physical ideas can be best illus- 
trated by contrasting the single and double sheet com- 
pounds. 

II. MODEL AND CALCULATIONS 

The simple model which contains the physics of bilayer 
systems is given by 

// = ^e„(k)cl(kK(k) (1) 

k,a 

+ \ ia,/3(k)cl(k)c^(k)-f /i.e., 

where k = {k^^ ky) is the in-plane momentum, a, f3 — 1,2 
represent the given plane indices and the sum over spins 
is implicitly included, i^./j represents the interplane cou- 
pling, through which electrons can hop from one plane to 
the other directly or through an intermediate state such 
as the chains or Y atoms. eQ.(k) are the band energies 
for the separate planes in the absence of an interlayer 
coupling. 

Due to the mirror plane symmetry of the bilayer, the 
above Hamiltonian can be diagonalized by considering 
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even and odd combinations of electron operators on the 
two planes. We introduce the new operators 

4(k) = i=[cl(k)±4(k)] (2) 

allowing us to rewrite the Hamiltonian in diagonal form 
H ^ ^[6+(k)c;(k)c+(k) +6_(k)cL(k)c_(k)], (3) 

k 

with the bonding and anti-bonding band energies given 
by e±(k) — e(k) ± tj_(k), respectively, where for sim- 
plicity we have taken degenerate bands in the absence of 
interlayer couphng: ei(k) = e2(k) = e(k). Also we have 
defined t±{k) = ii,2(k) = t2,i(k)- 

The intensity of Raman scatt^cing can be given as a 
scattering off an effective chargeE£l 

/5= 5] 7a,Mk)ci(k)c^(k), (4) 

where a,/3 are the plane indices and 7(k) is the Raman 
scattering amplitude. The effective charge is a noncon- 
serving quantity in contrast to the real charge. Here the 
Raman vertices 7 are related to the incident and scat- 
tering photon polarization vectors e^''^, resulting from a 
coupling of both the charge current and the charge den- 
sity to the vector potential. In general the Raman vertex 
depends non-trivially on both the incident and scattered 
photon frequencies. In the case of a bilayer system the 
light can interact with the electrons on the planes or the 
atoms in between. The effective charge can be rewritten 
in terms of the bonding and anti-bonding combinations 
of the bands, 

p = J2Uk)[clik)c+ik) + cUk)c_(k)] (5) 

k 

+7i,2(k)[c;(k)c+(k)-cUk)c_(k)]|, 

where 

7(k) =7i,i(k) =72.2(k) = i[7+(k)+7_(k)], (6) 
7i.2(k) = i[7+(k)-7-(k)]. 

This corresponds to a diagonal contribution which al- 
lows light to be scattered by density-like fluctuations on 
either plane 1 (71,1) or plane 2 (72,2) and an off-diagonal 
term (71,2) which allows light scattering on both planes 
simultaneously (see Fig. 2). Since in the limit of q — > 
the long-range intercell fluctuations are screened, only 
the intracell fluctuations remain. Therefore only charge 
transfer fluctuations between the atoms in the plane and 
outside the plane (e.g., Y and the chains), and between 
the planes can effectively cause light scattering. Since ex- 
perimentally the interlayer coupling is small, the result- 
ing vertex 71^2 must be smaller than the + combination, 



labelled as 7. Similar considerations have been applied 
to the case of multi-valley scattering in the conduction 
band of semiconductors. In the case when the intervalley 
scattering does not take place, no light scattering can oc- 
cur as the density fluctuations are screened completely. 
Scattering is restored when inequivalent valleys are con- 
sidered or the scattering is between different parts of the 
valence band. |— . 

According to Abrikosov and GenkinlHl if the energy of 
the incident and scattered frequencies uji,ujs are negli- 
gible compared to the relevant energy scale of the band 
structure, 7± can be expressed in terms of the curvature 
of the bands and the incident and scattered photon po- 
larization vectors e^'"^ as 

hm j^(k) = J2<?r^''^ (7) 

where terms of the order of 1 — ug/uji are dropped. This 
expression is valid given that the incoming laser light 
cannot excite band-band transitions. This effective mass 
approximation has been extensively used to calculate the 
Raman spectra in unconventional superconductors. 

However one might question the appropriateness of this 
approximation for the cuprates given that typical incom- 
ing laser frequencies are on the order of 2 eV - certainly 
on the order of the relevant electronic energy scale for a 
single band. For the case of two bands near or crossing 
the Fermi level, the effective mass approximation is even 
more questionable. Since in the bi-layer superconductors 
Bi 2:2:1:2 and Y 1:2:3 the band splitting of the even and 
odd bands is on the order of tens of meVs, the approxima- 
tion misses large terms corresponding to interband tran- 
sitions at relatively small energies. Detailed knowledge of 
the magnitude of the scattering amplitude in this case re- 
quires knowledge of the wave functions of the two states, 
which has currently not been investigated. Thereby, use 
of the effective mass approximation is uncontrolled in 
many band systems and comparisons to experiment in- 
tensities must be viewed with caution. A discussion of 
this point is elaborated in the Appendix. 

However, even if the above assumption does not hold in 
the actual experiments, this gives an insight into the rela- 
tive orders of 711 = 722 and 71^2 = 72,1, or in other words 
measures the splitting between 7+ and 7_ (see Eq. (6)). 
Without interlayer coupling the degenerate bands have 
the same dispersion such that 711 (k) = 722 (k), which 
holds generally due to reflection symmetry. 

An alternative approach is based on the experimen- 
tal observation that the spectra near optimal doping for 
a wide range of cuprate materials depends only mildly 
on the incoming laser frequency. Since the polarization 
orientations transform as various elements of the point 
group of the crystal, one can use symmetry to classify 
the scattering amplitude, viz., 

7(k; w/, W5) = ^7l(w/, W5)$L(k), (8) 

L 
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where $L(k) are either Brillouin zone (B.Z.H., orthog- 
onal over the entire Brillouin zone) or Fermi surface 
(F.S.H., orthogonal on the Fermi surface only) harmon- 
ics which transform according to point group transfor- 
mations of the crystalO. Representing the magnitude 
but not the k-dependence of both intra- and interband 
scattering, the prefactors can be approximated to be fre- 
quency independent and taken as model constants to fit 
absolute intensities. Thus we have simplified the many- 
band problem in terms of symmetry components which 
can be related to charge degrees of freedom on portions 
of the Fermi sheets. While sacrificing information per- 
taining to overall intensities, we have gained the ability 
to probe and compare excitations on different regions of 
the Fermi surface based solely on symmetry classifica- 
tions. This can be illustrated by considering the various 
experimentally accessible polarization orientations. 

Using an y coordinate system locked to the Cu02 
planes, incident and scattered light polarizations aligned 
along x + y,x — y for example transform according to Big 
symmetry, and thus 



$B (k) = cos(fc^o) - cos(fcya) -|- 



(9) 



where . . . are higher order B.Z.H. Likewise, e-^''^ aligned 
along x,y transforms as B2g'. 



*B2g(k) = sm{kxa) sin(fcya) 
The Aig basis function is 



(10) 



$^j^(k) = ao + a2[cos(A:a;a) -I- cos(A;j,a)] (11) 
+04 cos(fcj;a) cos(fcj^a) -I- oe [cos(2A;^a) + cos(2A;j^a)] -I- . . . , 

where the expansion parameters_ai determined via a fit- 
ting procedure with experimentll3. The Aig response is 
not directly accessible from experiments and must be 
obtained by subtracting several combinations of the re- 
sponse for various polarization orientations. 

By considering the k— dependence of the basis func- 
tions, it is clear that the Big part of the spectra essen- 
tially probes light scattering events along the or ky 
axes, i?2g probes the diagonals, and Aig is more of an 
average over the entire Brillouin zone. It is in this man- 
ner that information about the momentum dependence of 
the superconducting energy gap can be usefully extracted 
from the datatJ. That information lies in two important 
aspects: the low frequency power- law behavior of the 
spectra and the positions of the low energy peaks. 

At this point we can derive an expression for the cross 
section of the two band system. First we. define the 
k— dependent Tsuneto function A(k, lo;) aallu 



A(k, iu) 



A(k)2 



(12) 



X tanh 



2T 



2E{k) 



2E{k) 



Here ^^(k)^ = e(k)2 -I- A(k)2 and we have set ks = h = 1. 
We have neglected vertex corrections of the pairing inter- 
actions which while responsible for maintaining gauge in- 
variance and producing collective modestj, can be shown 
to have only a limited effect on the spectra at low energies 
for d— wave superconductorso. 

The Raman cross section is related to the Raman den- 
sity correlation function via the fluctuation-dissipation 
theorem. 



with 



l/T 



dre 



'{TApir)pl 



(13) 



(14) 



with Tr the time-ordering operator and the imaginary 
part is obtained by analytic continuation, icu uj + iO. 
In order to take into account the long-range Coulomb 
screening in the limit of g — > 0, the following formula 
must be calculated 



where x is given by Eq. (14), and 

»1/T 



(15) 



dTe-''^^Tr[p{T)p])=XlM^^)■, (16) 

l/T 



dre 



'{TAp{r)p]). 



(17) 







It is straightforward now to calculate the light scat- 
tering cross section for the two band model using Eqs. 
(5-6) and (12-17). We arrive at the following compact 
expression for the cross section: 



X 



;',(c.)=EWk)±7i,2(k)]2A'i(k,t 



(18) 



k,± 



Ek.±[7(k)±7i,2(k)]A±(k,c.) 
Ek,±-^±(k,^^) 



where A± is the Tsuneto function for the bonding/anti- 
bonding band, with energy gap A-i-(k). Rearranging 
terms, we can cast the result in terms of an addition 
of the result for each single band plus a mixing term: 

X'4(a;) - x+(^) + X" (^) + Ax"i^), (19) 

where 

XU^) = (20) 



E7l(k)A'i(k, 



Ek7±(k)A±(k, 
EkA±(k,^) 



4 



Ek7+(k)A+(k,c.) Ek7-(k)A_(k,c.)- 



(21) 



EkA+(k, 



EkA-(k,c.) 



Eqs. (19-21) are our central results. It shows that 
the total response can be considered as a sum of the 
single band contributions in addition to a mixing term 
which corresponds to odd combinations of fluctuations 
on the bands simultaneously (see Fig. 1). In order to 
discern the features of the mixing term, it is useful to 
write 7±(k) = 7° + A7±(k), where the average of A7(k) 
around the Fermi surface vanishes. The Fermi surface 
average of a quantity A(k) is defined as 



(^(k))F.s. 



Sk'^«k)-M)A(k) 



(22) 



Then it useful to consider the following cases: 

(1) If the layers are uncoupled, then 7+(k) ~ 7-(k) 
and the mixing term vanishes: Ax"{io) = 0. 

(2) If the scattering in the Aig symmetry channel 7 is 
independent of k for each band (as in the case of scat- 
tering on real charge), then Ax ~ (7+ — 7° )^ ~ 71^2 ■ 
This means that the light scattering induces an interlayer 
charge transfer. 

(3) Considering when the light induces charge transfers 
on different parts of the Fermi surfaces, A7±(k) must be 
included and includes different symmetry channels de- 
pending on the orientation of the incident and scattered 
light polarizations. Aig symmetry reflects the approxi- 
mate tetragonal symmetry of the Cu02 planes, while the 
part of 7 showing symmetries different than Aig (e.g.. 
Big, B2g and Eg) changes sign when symmetry transfor- 
mations are applied. The Tsuneto function depends only 
on the energy gap squared, and thus all the sums in Eq. 
(21) are zero except for the Aig component of^±{\i). For 
tetragonal materials, this means that the Big,B2g, and 
Eg channels do not have a mixing term and therefore 
the total response for these channels is simply the ad- 
ditive contributions from each single band. The mixing 
term only contributes to the Aig channel, the strength of 
which is roughly determined by the interlayer coupling. 

In the following section, we will investigate Eqs. (19- 
21) in detail for the single layer compound with different 
Fermi surfaces (section IIIA.) and bi-layer Y 1:2:3 (sec- 
tion IIIB). 



III. EVALUATION FOR MODEL 
SUPERCONDUCTORS AND DISCUSSION OF 
RESULTS 

We start by considering previous approachesi~0. The 
majority of the calculations have been performed by re- 
stricting the k-sum onto the Fermi surfaceH'uij. While 
giving qualitatively insightful description of the result- 
ing Raman spectra, this approximation does not capture 



important physics of charge dynamics over the entire 
Brillouin zone, such as the role of van Hove singulari- 
ties, and since the energy gap is a sizeable portion of 
the bandwidth, neglecting light scattering which breaks 
Cooper pairs at energies off the Fermi surface will miss 
large contributions to the overall cross section. This has 
been pointed out in the recent calculations of Branch and 
Carbottea. Moreover, since the Aig channel measures 
weighted averages over the Fermi surface, restricting the 
calculation for the Aig spectra to only the Fermi surface 
would generally overemphasis the effect of small changes 
of the Fermi surface manifold due to the presence of the 
gap singularities. Therefore, in our opinion the entire 
k-sum must be performed. 

In the following we confine ourselves to tetragonal sys- 
tems with in-plane lattice constant a which can be mod- 
eled in the absence of interlayer coupling by the following 
band structure: 

e(k) — ~2t[cos{kxa) + cos(fcj^a)] -I- 4i'cos(fc^a) cos{kya) 
-2t"[cos(2k^a) + cos{2kya)] - fi. (23) 

This form of the band structure results from a down- 
folding of an 8 baud Hamiltonian into a single band as 
described in Rcf.ta. The parameters tJl, and t" can 
either be taken from LDA calculationsEj or be taken 
as parameters—to fit the Fermi surfaces observed via 
photoemissionllSEl. 



A. Single layer case t± = 

In this section the T << Tc Raman response for a 
single CuO layer is obtained by evaluating the k-sums 
directly for the parameters for La 2:l:4o and Y l:2:3ll3. 

The Big and i?2g responses are plotted in Fig. 3 using 
an energy gap of d2;2_j,2 symmetry, A(k) — Ao[cos{kxa)— 
cos(fcj,o)]/2, with Aq chosen to be 12 (30) meV for La 
2:1:4 (Y 1:2:3), respectively, in order to matclt,the posi- 
tion of the Big peak to that seen in experiments. In both 
cases, the Big response for low frequencies varies as uj^ 
while the B2g varies linearly with uj. This can be shown 
to be generally true for these channels for any number of 
higher harmonics used, and results from a consideration 
of the density of states (DOS), which varies linearly with 
energy due to the nodes, and the behavior of the Raman 
vertex 7 near the gap nodes. Since the B2g vertex is finite 
near the nodes, the DOS determines the low frequency 
behavior. However, the gap and the vertex vanish at the 
same place for the Big channel, which in turn leads to 
the cubic behavior. This delicate interplay of vertex and 
energy gap can thus provide a unique determination of 
the nodal gap behavior. 

Moreover, additional information lies in the peaks of 
the spectra. A smooth peak is seen for B2g and a double- 
peak is obtained for Big. The lower peak is due to the 
energy gap while the other comes from the van Hove 
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singularity. The peak of the spectra due to supercon- 
ductivity in both cases is higher for the Big channel 
(wpeafe = '2Amax) than for B2g, (ujpeak ^ l.lAmax) since 
the gap maxima are located at the same part of the Fermi 
surface where the Big basis function is largest, and where 
the i?2g basis function vanishes. This is generally true 
regardless of Fermi surface shape, although the relative 
positions of the peaks can be affected slightly by chang- 
ing the^ underlying manifold (see Figs. 3a and 3b and 
RefsE^u). The shape of the spectra can also be slightly 
modified by considering higher harmonics of the-, energy 
gapB and/or by including final-state interactionsEJ or im- 
purity scatteringlJ. Apart from the last case, the power- 
law behaviors at low frequencies are unaffected and there- 
fore their observation is a robust check on the symmetry 
of the energy gap. 

The van Hove peak only shows up in the Big channel 
since the B2g channel assigns no weight to the location 
of the van Hove points. The van Hove peak lies at a 
higher energy determined by its location off the Fermi 
surface and the value of the energy gap there. There are 
indications that a double-peak featuije has been seen in 
underdoped Y 1:2:3 with a Tc = BOifQ. As the van Hove 
moves further away from the Fermi surface for example 
with doping, the peak shifts to higher frequencies and has 
a smaller residue. Inelastic scattering at higher energies 
will also act to smear out this feature and thus a single 
smeared and perhaps asymmetric peak would be seen if 
the van Hove lies quite near the Fermi surface, and could 
in principle move the Big peak slightly above JiAj^ax- 
This may be reflected in the experiments of Ref.Q which 
saw a sensitivity of the Big peak position to doping. 

In Fig. 4 we show a calculation of the unscreened 
and screened Aig response for a d^^-y^ paired super- 
conductor for the Y 1:2:3 parameters. We have cho- 
sen the coefficients in Eq. (11) using the effective 
mass approximationllj. By again considering the nature 
of the nodal structure of the energy gap, since the Aig 
response measures a weighted average of charge fluctu- 
ations around the Fermi surface, it can be shown that 
the low energy part of the Aig response must vary with 
energy in the same way as the DOS. 

The screened and unscreened responses at higher ener- 
gies are very different. For the unscreened response, two 
peaks arise corresponding to the singular contribution of 
pair breaking along the axes and the van Hove contribu- 
tion as in the Big channel. Screening completely reorga- 
nizes the spectra to remove both peak contributions, re- 
placing them with a much smoother function with more 
spectral weight at lower energies and a much reduced 
spectral weight at higher energies. This is due to the 
screening function (2nd part of Eq. (15)) which also con- 
tains peaks at 2Amax due to the weighting along the 
and ky axes, and the van Hove. These peaks cancel the 
unscreened peaks and create a mild shoulder lying near 
2Ajnax, as seen in Fig. 4. 

These results are different than those obtained previ- 
ously for the case of a simpler evaluation of Eq. (18) 



on a cylindrical Fermi surfaceEl and for a more thorough 
evaluation of Eq. (18) using the full k-sum and a consis- 
tent treatment of the harmonics for both the band struc- 
ture, Raman vertices, and energy gapcl. Since publica- 
tion of Ref.ld, a numerical error was detected in the code 
which once corrected also produced curves similar to the 
one shown in Fig. 4. Thia.is treated and corrected ip. 
a forthcoming publicationEEl. We also note that Ref.cl 
also produced an Aig response which is peaked near 2A, 
but remark that these calculations were performed by 
neglecting the real parts of Eq. (20). 

In an effort to understand these differences we have in- 
vestigated Eq. (20) for various choices of band structure, 
Raman vertices, and additional harmonics added to the 
energy gap. While we found that the linear power-law 
rise of the Aig response is robust to changes in parame- 
ters, the position of the peak in the spectra can vary any- 
where less than or equal to 2A depending on the choice 
of parameter sets. This suggests that the calculations for 
the Aig screened response are unstable to the inclusion 
of higher order Brillouin zone harmonics for either the 
energy gap, band structure, or the Raman vertex. 

Specifically, we have included the next few harmonics 
for a d^2_y2 energy gap, 

A(k) = Ao| [cos(fca;a) - cos(fcya)]/2 + (24) 
^[cos(fcj,a) - cos(fcj^a)]^ -|- ^[cos(fcj,a) - cos(fcj,a)]'''|, 

and used the parameters ai iA-§q- (H) for the vertex de- 
rived from the band structurellj. The results are summa- 
rized in Fig. 5, which shows the various forms for the Aig 
screened response using different combinations of energy 
gap harmonics. This result suggest that the weighting 
of the Aig vertex in the Brillouin zone depends strongly 
on a variety of parameters and can easily weight regions 
where the energy gap is quite different for small changes 
in the parameters. This results in varying positions of 
the peak. In general, including higher harmonics for the 
energy gap pushes the peak in the Aig channel downward 
from near 2Amaxi which is in agreement with the results 
obtained by Branch and Carbotteo. We are then to con- 
clude that the position of the Aig peak can not provide 
useful information without detailed knowledge of the fine 
structure of the energy gap itself even for the case of a 
single band near the Fermi level. 

In this manner, similar symmetry considerations can 
be made for various types of energy gaps. The Raman 
spectra calculated for various energy gaps is summarized 
in Table 1 and compared to the low energy DOS N{uj). 
The dominant contribution to the Raman lineshape is 
due to the location and behavior of the energy gap near 
both the nodes and the maxima. Subsequently, different 
energy gaps produce different line shapes as summarized 
in Table 1. Here it is important to note that the cu- 
bic rise of the spectra in any channel requires that both 
the vertex and the energy gap have the same behavior 
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near the nodes. In tetragonal systems, this requires the 
presence of a dx2-y2 energy gap. By considering small 
orthorhombic distortions, the Aig and Big channels be- 
come mixed and therefore a linear rise with frequency 
could be obtained at low energies in a region determined 
by the amount of symmetry breaking. Since experimen- 
tally the low frequency part of the Big spectra rises cu- 
bically in the most tetragonal systems, while in Y 1:2:3 
the spectra may have an intrinsic small linear part, this 
strongly suggests the likelihood of a dx2-y2 pair state 
for these systems. Exploration of the low frequency part 
of the spectra could put stringent constraints on other 
pair-state candidates. 

B. Bi-Layer 

In this section we consider the bi-layer superconductpi: 
Y 1:2:3. For the interlayer couphng t±, we follow Ref.E£l 
and define 

t^Ck) = t^[cos{kxa) — cos(fcya)]^, (25) 

with t±/t — 0.2. In the absence of inter- layer coupling 
the chemical potential is chosen so that (n) = 0.8. The 
two Fermi surface sheets are shown in Fig. 6 for these 
choices of the parameters. Lastly, we choose the energy 
gaps on the two bands to be of the dx2-y2 type, A±(k) — 
A±[cos{kxa) — cos{kya)]/2, and further choose A± = 30 
meV for both bands for simplicity. 

The resulting spectra for the bonding and anti-bonding 
bands, as well as the mixing term for Aig are shown in 
Fig. 7. Here again we have used the effective mass ap- 
proximatiotufor the expansion parameters of Eq. (11) for 
illustrationEj, although as indicated earlier this approxi- 
mation does not hold for this bi-layer system. Consider- 
ing the Big and B2g channels, which are simply additive, 
it can be inferred from experimental result which shows 
only one peak in these channels that either the energy gap 
is nearly identical both in symmetry and magnitude or 
alternatively that the Raman response is given predomi- 
nantly by one band. Otherwise two distinct peaks would 
appear (if gaps have different symmetry and/or magni- 
tude) and identical power-law behavior for each channel 
would be seen (if the_jgaps were of different symmetry 
but same magnitude)E3. Since the band splitting is rela- 
tively small in Y 1:2:3 and even smaller in Hi 2:2:1:2, the 
second possibility in our opinion is unlikelyO. Therefore, 
the energy gap on each band must be similar on energy 
scales determined by the amount of inelastic scattering 
which would smear a double-peak feature. Again we note 
that the sharp nature of the Big peak would imply that 
its position and shape would be most sensitive to the 
location of the van Hove points fQc-|each band and the 
maximal value of the energy gapaB'El. In fact, as shown 
in the Fig. 8, the van Hove peak in some cases can be 
stronger than the 2Amax peak which could erroneously 
lead to an overestimated value of 2Amax/T^c- 



Due to the small bi-layer splitting the mixing term for 
Aig is small, as shown in Fig. 7. Two very small peaks 
with opposite curvature at roughly 2Amax for each of 
the bands are seen due to the odd combination of inter- 
layer charge transfer, as suggested by Krantz and Car- 
dona. These peaks are non-divergent and are strongly 
suppressed compared to the contributions from the in- 
dividual bands. Since adding electron coHrfijlations sup- 
presses the bi-layer splitting even furtherEja the mixing 
term is only of minor importance and the resulting spec- 
tra can be well approximated as the sum of the contri- 
butions from the two individual bands. In this way, a 
consistent picture can emerge since experimentally the 
electronic Raman spectra do not djifer substantially in 
the bi-layer or single-layer materialsoO. 

We have also checked the effect of abandoning the effec- 
tive mass approximation for the vertex as well as adding 
additional harmonics to the energy gap, as considered 
in the previous section. Our results are shown in Fig. 
8 for the parameters chosen to most closely model the 
experimental dataO: 

Big : 7(k) ^ cos{kxa) — cos(fcya), 
B2g : 7(k) ~ sin(fc3;a) sin(fcj^a) — [sin{kxa) sin(fcj,a)]'^, 
Aig : 7(k) ~ a4 cos(fca;a) cos(fcj,a)-|- 
a6[cos(2fca;a) + cos(2fcj^a)]; 04 = —2t' — t^.a^ = At" -\- 1±, 

A2 = 0.5, Ao = 30meV. 

The results show that the peak of the Aig spectra moves 
downward in frequency with the inclusion of the addi- 
tional energy gap harmonic as for the single band case, 
while the mixing term and power-law behavior remains 
unaffected. The B2g peak is pushed to lower frequencies 
by the inclusioH of the higher harmonic for the energy gap 
as seen in Ref.EI, but can be shifted back by including the 
next harmonic for the Raman vertex. The resulting peak 
lies roughly 30 percent lower than the peak in the Big 
spectrum. The van Hove feature does not change appre- 
ciably and thus the Big channel changes are quite small. 
The power-law behavior is unaffected, as remarked ear- 
lier. From these results we can then conclude that the 
Raman response for the bi-layer can be well modelled by 
the added response calculated for the two separate en- 
ergy bands with the only difference is that the use of the 
effective mass approximation is invalid (see the Appendix 
for more details). 

IV. CONCLUSIONS 

We have seen that electronic Raman scattering con- 
tains a wealth of information concerning the nature 
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of superconductivity in unconventional superconductors 
which can be extracted from symmetry considerations 
alone. When applied to the experimental data for 
optimally doped cuprate materials, good agreement can 
be obtained for dx2_y2 pairing in single band materials, 
and nearly identical gaps of d^i^yi symmetry for bi-layer 
systems. The analysis is based on the prominent behav- 
ior of the spectra at low frequencies and the approxi- 
mate position of the maxima for each channel. It is also 
worthwhile to point out that similar considerations can 
be applied to the optically active phonon lineshapes be- 
low Tckll. 

In this paper we confined ourselves to the simple case of 
clean, weakly coupled unconventional superconductivity. 
As pointed out earlier, like other correlation functions 
Raman in clean systems is insensitive to the phase of the 
order parameter around the Fermi surface. Moreover, 
as a consequence of taking the g — > limit, the consid- 
ered spectra vanishes as one approaches T^, due to phase 
space constraints. Both of these deficiencies have been 
remedied receiitJ5i-.by considering both the effect of impu- 
rity scatteringlH^cd. As shown previously, disorder effects 
on unconventional as opposed to conventional supercon- 
ductors can in principle be used to determine whether 
the gap averages to zero around the Fermi surfaceE3. 
In particular, the disorder dependence of the exponent 
of the low frequency and/or temperature dependence of 
the Big chajMijel provided a useful signature of gap phase 
propertiesEIlE3. While the resulting spectra could be ap- 
plied to low frequencies, the flat "Marginal" behavior at 
high frequencies cannot be reproduced by this method. 
However, that can be corrected by considering spin fluc- 
tuatioru.4(:a,tl£jring or electron-electron scattering due to 



nestingllljijO. However, more work is needed to cor- 
rectly describe the full channel dependent Raman spectra 
for large frequency ranges for optimally doped cuprates. 

Further, we have seen that the lineshape of the Aig 
screened response is sensitive to parameter changes and 
to the number of harmonics used for the energy gap. 
However, the low frequency behavior and the general 
"flatness" of the lineshape seems to be independent of pa- 
rameter choices. We remark that a consistent treatment 
of the number of harmonics used for the band structure, 
Raman vertices, and energy gap must be considered to 
give meaningful results for the peak position. Without 
a more detailed knowledge of the actual k-dependence of 
the energy gap (symmetry is not enough), a theoretical 
description of the peak of the spectrum in the Aig chan- 
nel is incomplete. In addition, a recent study found that 
the position of the Aig peak did not move uii|dcr pres- 
sure, in contrast to the peak in the Big channeled, which 
may indicate that the peak in Aig may not be directly 
related to the charge degrees of freedom. Clearly, further 
work is needed. 

Moreover, further work is needed in order to under- 
stand the prop|eijti£a of the Raman spectra away from 
optimal dopin^a~Q'E£l. It is suggested that the Big peak 
may be the most sensitive to doping due to (1) the sharp- 



ness of the peak, (2) the role of the van Hove singularity, 
and (3) the sensitivity of the Fermi surface along the 
kx , ky axes to small changes of the band structure pa- 
rameters. However, in order to understand the effects 
of doping and to arrive at a more complete picture of 
Raman scattering in the cuprates, a particular model of 
superconductivity is needed as is a more sophisticated 
approach for handling the spin degrees of freedom. 
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V. APPENDIX: EXTENSION TO BANDS 
CROSSING THE FERMI LEVEL 

In this appendix we develop formulas for electronic Ra- 
man scattering applicable for a multiband model. Our 
goal is to show that 

i) the use of the effective mass approximation for the 
Raman vertex is unjustified in case of more than one 
band crossing the Fermi level, and 

ii) there can be no divergence in the Aig Raman spec- 
trum in the superconducting state. 

First we recall that the differeniial electronic Raman 
scattering cross section is given byllj 



(Pa 
dflduj 



1 



-hui/keT 



ImXp+,p(w), (26) 



where a; = — is the difference between the incoming 
and scattered light frequency, rg = /mc^ is the Thom- 
son radius, and 



xp+.pW = r0W([/5+(q'*).^(q'O)])- 



(27) 



For a Bloch electron system with one-electron states 
|n,k) the effective density operator is represented as 



(28) 



n,n' k.cr 
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with q = — is the difference of the incoming and 
scattered hght wavenumber. In the experimental situa- 
tion relevant for Raman scattering (j*,(7^ ^ kp, and the 
Raman vertex simplifies tcO 



(29) 



m ' ^ 



(n,k|e''p|j/, k)(i/,k|e^p|n',k) 



e„'(k) -e^(k) + 
(n, k|e*p|i/, k) {y, kje-'plTi', k) 



e„'(k) - e,,(k) - hbj^ 



where e* and e* are the polarization vectors of the in- 
coming and scattered light, p = —iTiW, and £n(k) is the 
energy of the Bloch electron. The correlation function of 
the effective density is now readily evaluated as 

n,n' k,(7 

/[£„(k)]-/[£„,(k-q)] 
hcu + iS — en(k) + e„'(k — q) 

The above susceptibility consists of both intraband (n = 
n') and interband {n ^ n') contributions. 

In the context of an N band model it is implicitly 
assumed that out of the infinite bands of a Bloch sys- 
tem only N bands (with typical bandwidth W) cross the 
Fermi level, while all the other bands are far above or 
below those N bands of interest, separated by a gap G 
of order G ^ W. This is the largest energy scale of the 
problem, therefore when we consider Raman scattering in 
this model we must assume Tilu'^, Tilo^ <^ G, although the 
energies of the incoming and scattered light may well ex- 
ceed those of the optical transitions between bands close 
to the Fermi level (non resonant scattering) . Under these 
assumptions only terms with indices n and n' both label- 
ing bands crossing the Fermi level will contribute to the 
Raman susceptibility in Eq.(30). The relevant Raman 
vertex will further simplify to 



7„,„'(k) = e'e^5„^„' + 



(31) 



E 



(n,k|e"p|i/, k)(iy,k|e''p|n',k) 



(n, k|e*p|j/, k)(t/, kje^'pln', k) 



[£„.(k)-£,(k)]- 



where the prime on the summation indicates that v runs 
only over indices belonging to the faraway bands. Since n 
and n' belong to bands close to the Fermi level, the light 
energies are negligible compared to the large band gaps in 
the denominators of Eq. (30). On the other hand when v 
assumes values belonging to the " close" bands, the band 
gaps are negligible compared to the light energies, and 
it is easily seen, that those terms yield contributions of 
the order of WTiu) /fiuj^fLLd'^ « oj/uj^'^. Under usual ex- 
perimental circumstances {lo <C uj^) these terms can be 
neglected in the Raman vertex. 



Since we are interested in the low frequency response 
and interband transitions contribute to ImXp+ ^ only 
above some finite threshold, in the following we restrict 
ourselves to intraband {n — n') transitions. The Raman 
susceptibility is now given by 



Xp+,p(^) 



n,k,cr 



[74k)]= 



/[£„(k)]-/[£„(k^q)] 
hu + iS - £„(k) -I- £„(k - q) 



(32) 



where the Raman vertex is 



J 1 A (n,k|e*p|j/,k)(i/,k|e'p|n,k) 4- c.c. 
74k)=ee £„(k)-£.(k) 



(33) 

real, and an even function of k. The susceptibility is 
additive for those bands crossing the Fermi level, however 
the Raman vertex reduces to the familiar effective mass 
formulsO 



-7 fkl - Vc- m 5^£„(k) 

a, [3 



(34) 



only for the case of a single band crossing the Fcmii level. 
This is because in the effective mass formulaEZi v runs 
over all values except n, while for more than one bands 
crossing the Fermi level there are some other terms also 
excluded from the sum for the Raman vertex in Eq.(33). 
The missing terms are in general of order unity due to the 
relatively small gaps between bands close to the Fermi 
level. In the Raman vertex however, these terms are 
negligible because of the large laser light energies (see 
Eq.(29)). Therefore in a multiband problem knowledge 
of the electronic energies £„ (k) is not sufficient for a pre- 
cise evaluation of the Raman vertices 7„(k). Using the 
effective mass formula in this situation renders the anal- 
ysis qualitative at best. 

Consider our system in the superconducting state. As- 
suming that superconductivity develops independently in 
the bands crossing the Fermi level, the intraband Ra- 
man response will again be additive, and the well known 
formulali3 is generalized as 



Xp,p(^) 



n,k 



:7„(k)]2A„(k,c.), 



where 
A„(k, w) 



tanh 



£^n(k) 



2kBT 



4|A„(k)|Vg„(k) 
AEl{\<i) - {hoj + isy 



(35) 



(36) 



is the Tsuneto function for the n*'* band involving the su- 
perconducting order parameter A„(k) and quasiparticle 
energy £'^(k) = £2 (k) | A„(k)p of that band. Although 
it is not trivial, it can be shown that for the N band 
model outlined above the Raman vertex derived for the 
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normal state in Eq. (33) is an excellent approximation 
for the one to be used in the superconducting stateCJ. 

The screening effect of the long range Coulomb inter- 
action can be taken into account by including coupling 
of the effective density to the density {p) fluctuationsEa. 
In the long wavelength limit appropriatQ-£o^ Raman scat- 
tering the screened response is given bj 



scr 
Ap,p 



Xp,/ 



XP;PXP;P 

Xp,p 



with 



and 



XpA^) 



^7„(k)A„(k,w) 



Xp.pi^)^ 



(37) 



(38) 



(39) 



Note that there is no interband contribution to Xp-p ^-^d 
Xp,p in the long wavelength limit appropriate for the 
present analysis. On the other hand, the second term 
in Eq. (37) mixes the bands (unless Xp,p — 0)i thereby 
making the screened response non additive. The screen- 
ing effect of the long range Coulomb interaction can be 
made more transparent by rewriting Eq. (37) with the 
help of Eqs. (35), (38) and (39) as 



by screening, since 7 at w = 2|A|ma2; takes exactly the 
same value as 7(k) does at the critical k points. This 
is because for uj = 2|A|„ia2; the Tsuneto function assigns 
infinite weight in 7 to exactly that 7 value. Therefore 
again, any peak in the experimental Aig spectrum is not 
directly related to the superconducting gap. 

In a multiband situation the question might arise as to 
whether the above one band result holds for the removal 
of any divergences in the Aig symmetry. Fueling doubts 
about it is the fact that in a one band model a momentum 
independent Raman vertex would completely be screened 
out leading to no response at all, while in the multiband 
case different but still momentum independent 7„ values 
will obviously lead to finite response according to Eq. 
(40). The existence of a divergence however is a more 
subtle problem. Any divergence can only come from one 
of the Tsuneto functions at one of the different frequen- 
There is no reason to assume that 



2|A 



n \ max ■ 



xfji^) 



[7„(k)-7(c.)]2A„(k,(^), 



(40) 



the maximum gap is the same for any two bands if those 
bands are different. For a given a;„ only A„ can give rise 
to divergence, the contribution of all the other bands in 
Eq. (40) is finite at that frequency. Similarly, j{wn) will 
be determined solely by the divergent contributions from 
the n*'' band both in the numerator and denominator of 
Eq.(41). Therefore as far as the existence of divergences 
in the Aig spectrum is concerned, the multiband problem 
reduces to the one band problem with the same result: 
any possible divergence at the various gap maxima in the 
Raman spectrum is removed by screening. 



where the average of the Raman vertex (which is screened 
out) involves all the relevant bands and depends on fre- 
quency as 



En,k7»(k)A„(k,Lj) 

E,, k •^n(k,w) 



(41) 



The Raman response as given by Eqs. (40) and (41) 
has been thoroughly investigatecLfor the case of only 
one band crossing the Fermi levelEJ. It was shown, that 
for tetragonal symmetry and for a quasi two-dimensional 
electron system only the Aig response is screened. More- 
over, since the Tsuneto function picks up divergent con- 
tributions only for uj = 2|A|maa;; and those contributions 
come from the neighborhood of points on the Fermi sur- 
face where the maximum value of | A| is attained, the nec- 
essary condition for divergent Raman response is that the 
Raman vertex does not vanish at those points. For exam- 
ple, in case of a superconducting gap of dx'^-yi symmetry, 
the Big response is divergent at w = 2|A|ma£c translat- 
ing into a robust peak at that frequency in experiments, 
while there is no divergence (anywhere) in the i?2g re- 
sponse. Therefore the frequency of any observed max- 
imum in that spectrum is not related to the supercon- 
ducting gap maximum. The unscreened Aig response is 
also divergent at 2|A|mQa;, but this divergence is removed 
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TABLE I 



Summary of Raman response for various pair state candidates for clean, tetragonal superconductors.* 
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* Here phase is undetermined and thus | dx^-y2 \ would yield the same as dx2_y2. Moreover, orthorhombic distortions 
will mix Big and Aig channels. 

Parameters chosen to split node at 45° into two modes at 45° ± 5° and a subsidiary maximum equal to 20 percent 
of gap maxima. 

Parameters chosen such that (n) = 0.8. 

Alg peak position depends sensitively on harmonic representation of band structure, vertices, and energy gap (see 
text for a discussion). 
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FIG. 3. Big and B29 Raman responses for a single band plotted for parameters sets for (a) La 2:1:4 and (b) Y 1:2:3. Here a 
filling (n) = 0.8 and t = 100 meV, and A(A;) = Ao[cos(fca;a) — cos(A:j,a)]/2 are used for both figures. The other parameters are: 
(a) t'/t = 0.16, t" = 0, Ao = 12meV; (b) t' /t = 0.2,t"/t = 0.25, Ao = 30 meV. The small peak at higher frequencies in the Big 
channel is due to the van Hove singularity. 
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FIG. 4. Aig response calculated for a single band with and without Coulomb screening for the Y 1:2:3 parameters. 
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FIG. 5. Comparison of the screened Aig response obtained for different number of gap harmonics used than for Fig. 4 (which 
is redrawn here and labelled as (a)). All parameters used are the same as Fig. 4 except the following: (a) no changes, (b) 
Ao = 32 meV, Ai = 1, (c) Ao = 35 meV, Aa = 1 (d) Ao = 36 meV, Ai = Aa = 1. 
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FIG. 6. Fermi surface for the absence (labelled as 0) and presence of interlayer coupling (labeled ± for the bond- 
ing/anti-bonding bands). 
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FIG. 7. Raman response calculated for Y 1:2:3 bi-layer for various chaiuicls as indicated. The +( — ) indicates the bonding 
(anti-bonding) band, respectively, the symbol M is the mixing term [Eq. (21) which only contributes for the Aig channel], and 
the total response is indicated by the symbol T. The peaks are located a different positions for the Big channel due to the 
different energy dispersions of the two bands. The van Hove peak lies at much higher energies for the bonding band due to its 
larger separation from the Fermi level. 
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FIG. 8. Raman response calculated for Y 1:2:3 bi-layer for various channels as indicated for different parameter sets than 
those considered in Fig. 7. See text for details. 
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